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We explore a dark-matter model in which there are two dark-matter species nearly degenerate in 
mass, with e = AM/M <C 1. The heavier particle undergoes two-body decay with a half-life T, 
to the lighter dark-matter particle and a noninteracting massless particle. Unlike previous work 
on decaying dark matter, we explore the regime x > 100 Myr and non-relativistic kick speeds 
vjt/c — e. Using a set of N-body simulations of isolated dark-matter halos, we show how halos 
change as a function of X and Vfc. We find that X < 40 Gyr is ruled out for Vk > 20 km s _1 
(e > 10~ 4 ) when we compare the simulations to observations of dwarf-galaxy- to cluster-mass 
dark matter halos. We highlight which set of observations should provide better future constraints 
for decays and other types of dark-matter physics. 
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1. Introduction 

Most of the gravitationally attractive mass-energy of the Universe consists of dark matter, but 
the physics of dark matter is poorly characterized [jl|]. The most popular dark-matter candidate 
class, the weakly-interacting massive particle (WIMP), consists of particles that are heavy (m > 
100 GeV), self-annihilating, stable, and have low interaction probabilities with other particles [Q]. 
However, there are a number of other well-motivated candidates with different phenomenology 
[§, 0]. Dark matter may have strong elastic self-interactions, it may not self-annihilate, or may be 
unstable ^, R]. The identification of dark matter depends on uncovering its phenomenology. 

In this proceeding, we discuss constraints on the stability of dark matter. In particular, we 
consider a dark-matter model that consists of two nearly degenerate WIMP species, X and Y, in 
which X decays to Y and a (nearly) massless particle on cosmological time scales. This is one of 
the simplest extensions to the WIMP model. The mass difference e = AM /M is related to the recoil 
speed of the daughter Y particle v^/c = £ if £ -C 1. There are thus three parameters for this model: 
Mx, the mass of the X particle; Vk (or £), the recoil speed of the daughter Y; and T, the half life. 

While constraints on this model are stringent if the (nearly) massless particle in this model 
is a gamma ray or decays to light leptons there are few constraints if the (nearly) massless 
particle has negligible interactions or is a (or decays to) neutrino. To constrain a generic nearly- 
degenerate decay model, we focus on the structure and evolution of dark-matter halos. Decay can 
affect halos in two ways: daughter particles can escape, causing mass loss; and the dark-matter 
density decreases within the halo in response to mass loss or, in the case the daughter particle 
remains bound, to the difference in kinetic energy between the parent and daughter particles. In 
practice, this means that the value of Mx is irrelevant as long as it is similar to WIMP-range masses. 

We use simulations to characterize the evolution of dark-matter halos in response to decays, 
as a function of and t; these are described briefly in Sec. || and in more detail in Refs. H |10| ]. 
In Sec. |3| we use the simulations as a map the decay parameters to observables in order to use 
the properties of dark-matter halos inferred from observations of galaxies, groups, and clusters to 
constrain the decay parameter space. We discuss our key findings in Sec. |j. 

2. Simulations 

We simulate isolated equilibrium dark-matter halos using the parallel A-body code GADGET- 
2, modified by us to handle decays JTT| , ^]. Since we consider % to be long, we assume that halos 
collapse before there is significant decay. Thus, the initial conditions for the halos will be well- 
described by the expected properties of cold-dark-matter (CDM) halos. Our halos are spherically 
symmetric, and the density p(r) as a function of galactocentric radius r is taken to be the Navarro- 
Frenk- White (NFW) form [||], 




(2.1) 



where p s is the scale density, and r s is the scale radius, such that the halo concentration is 
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where is the virial radius of the halo. We define virial quantities with respect to the spherical 



top-hat overdensity, as described by Ref. [|13fl. 



The initial velocity distribution must be chosen carefully in order for the halo to be in dy- 
namical equilibrium in the absence of decays. We take the velocity ellipsoid to be isotropic, a 



reasonable first approximation for dark-matter halos []14|]. The procedure to initialize the velocities 
is described in more detail in Ref. [^]. 

We begin every simulation with identical initial conditions, with 10 6 particles in the virial 
radius, varying only the concentration c of the NFW profile, particle lifetimes, and kick velocities. 
We simulate virial mass M v ; r = 1O 12 M halos only, but we can extrapolate our results to other halo 
masses by considering the structural changes to the halos as a function of the decay parameters 
with respect to the halo virial speed and dynamical time. We perform a set of 100 simulations, with 
v k = 10, 50, 100, 200, and 500 km s" 1 and % = 0.1, 1, 10, 50, and 100 Gyr, and with c = 5, 10,20 
and 30. Since the virial speed of the halo is v v i r = y/GM~jR~r ss 130 km s" 1 using a spherical 
top-hat overdensity definition of virial parameters, the kick speeds were chosen to bracket the virial 
value. The simulations are discussed in depth in Refs. [T(J . 

3. Results 

In this section, we show what kind of constraints one can set on the decay parameter space 
given the observations described below. Those interested in the details of the evolution of the 
dark-matter halo density and velocity-dispersion profiles as a function of decay parameters should 
consult Sec. Ill of Ref. 

3.1 Cluster Mass Function 

The comoving galaxy-cluster mass function dn/dM is a sharply falling function of halo mass 
at all redshifts, and is a probe of both the Hubble expansion (via the conversion of redshifts to 



comoving volumes) and the growth function of structure Q15Q. In the absence of decays or modifi- 
cations to gravity, the cosmological parameters inferred from probes of homogeneous cosmology 
(e.g., the Hubble expansion) should agree with those inferred from probes of the growth of struc- 
ture, and those parameter values should be consistent across cosmic time. However, decays change 
the growth of structure, with mass loss from halos becoming increasingly important with time. If 
the decays are non-relativistic, then the Hubble expansion should not deviate from what would be 
predicted from the cosmic microwave background, but the growth function would be different. 

To constrain decays using the cluster mass function, we take the highest-normalization z = 
cluster mass function allowed by high-z probes of cosmological parameters (e.g., the cosmic 
microwave background [[I]]) assuming a ACDM cosmology, and map the CDM halo masses to 
halo masses after a Hubble time of decay using our simulations. In order for a point in decay 
parameter space to be allowed, the z = cluster mass function from the decay simulations must 
lie above the minimum cluster mass function allowed by X-ray and optical observations of low- 
redshift clusters [16]. We show an example of constraints in the left-hand plot of Fig. [I], in which 
the comoving number of clusters above mass M vu is plotted against M v - ir . The colored lines indicate 
(top to bottom) the highest cluster mass function consistent with WMAP seven-year data, the mean 
cluster mass function predicted from WMAP, and the 1— a lowest mass function consistent with 
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Figure 1: Examples of constraints on decays from {left) the z — cluster mass function and (right) the z = 
mass-concentration relation. See text for details. 

z = observations of clusters Jlr^]. The dot-dashed lines indicate cluster mass functions for Vk = 
2000 km s _1 and (top to bottom) t = 100,40,20, and 5 Gyr. In general, T > 30 Gyr for v k > 
1000 km s _1 (£ > 10~ 2 ) in order to be consistent with the observed cluster mass function. 



One can also use the redshift-dependent cluster mass function to probe decays (see, e.g., [17]). 
The strength of this approach is that the deviations of the cluster mass function from a ACDM 
model chosen to match z = observations should increase with look-back time. The challenge of 
this approach is that the redshift-dependent cluster mass function will also probe non-A models for 
the accelerated expansion, and there may be degeneracies between decay and modified gravity in 
the redshift evolution of the cluster mass function. 

3.2 Mass-Concentration Relation 

The mass distribution within dark-matter halos depends on the formation time; small halos 
that form early when the Universe is more dense have a higher central-density-to-virial-density ratio 
than larger halos that formed later, causing a mass-concentration relation [ff3]. High-normalization, 
high-il m cosmologies tend to produce higher concentrations for fixed mass than for low-normal- 
ization, low-£2 m cosmologies because the typical formation time of halos is earlier. Decays lower 
the concentrations of dark-matter halos due to the negative heat capacity of self-gravitating systems, 
so the observed mass-concentration relation should provide good constraints on decay. 

We show how decays affect the mass-concentration relation in the right-hand plot of Fig. ||. 



On this plot, the best-fit NFW concentration (Eq. 2.2) is plotted against the virial mass of halos at 



z = 0, with halos spanning the range of those hosting L* galaxies to massive clusters. The upper 
and lower solid lines show the predicted mass-concentration relation for WMAP one- and three-year 
cosmologies, respectively, the former of which is approxately the 2— a upper limit from WMAP 
seven-year data [jl^]. The short-dashed and dotted lines show estimates of the mass-concentration 



relation from strong-lensing and X-ray observations of clusters [20]. The selection function of 
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z re = 7 satellite + large-halo limits 
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v k [km s" 



Figure 2: Exclusion limits in the z — Vk parameter space. See text for details. 

the observations is poorly understood but is generally thought to prefer high-concentration, early- 
forming relaxed clusters. The magenta-shaded region represents the 1— a uncertainty in the central 



mass-concentration relation from weak lensing assuming a ACDM cosmology Q21[]. The error bar 
in the top-left corner of the plot indicates the scatter in the simulated mass-concentration relation 
and the typical I— a uncertainty in the fit for the mass-concentration relation from observations. 

To find the most conservative constraints on decay, we assume a high-normalization, high-il m 
Universe consistent at the 2— a level with the WMAP seven-year analysis. We consider a point in 
decay parameter space to be ruled out if its corresponding mass-concentration relation lies below 
the magenta-shaded region of Fig. |L since weak lensing has a better-understood selection function 
than strong lensing or X-ray clusters and groups. In the right-hand panel of Fig. [|, we consider 
Vk = 500 km s _1 , with the black dashed lines corresponding to (top to bottom) x = 100,40,20, and 
5 Gyr. In general, we find x < 30 Gyr to be excluded for Vk > 200 km s _1 . This is more exclusive 
than for the cluster mass function because weak lensing probes much smaller dark-matter halos, 
and decays have the greatest effect on halo structure if Vk > v v - ir . 



3.3 Local Group Dwarf Galaxies 

The Local Group has traditionally been a source of inspiration for dark-matter theories, es- 
pecially since by 1999, simulations with ACDM cosmologies indicated that there ought to be far 



more subhalos in the Milky Way halo than there were satellite galaxies observed [22]. However, the 



picture of Local Group satellite galaxies has changed drastically in the past few years. The number 
of known Milky Way satellites has more than doubled, and their kinematic properties have been 



carefully characterized Q23| , g4J. Critically, all known Milky Way satellites have high dark-matter 
density; the mass within 300 pc of the centers of these satellites is M300 >5 x 10 6 M Q regardless 
of luminosity. Any explanation for Milky Way satellite population must correctly reproduce the 
central density of the satellites as well as the luminosity function and the number of satellites. 

In Ref. JTfj|], we used a hybrid technique to construct realizations of Milky Way satellite pop- 
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ulations as a function of decay parameters. We used analytic merger trees for a high-normalization 
ACDM cosmology in order to find subhalo populations for Milky Way-mass halos. We used the 
simulations described in Sec. || to map the properties of the ACDM halos and subhalos to those at 
z = assuming specific Vk and z. Finally, we populated subhalos with stars according to a crude 
star-formation prescription tied to the redshift of reionization z re , for which we take a fiducial value 
z re = 7. This star-formation model is meant to illustrate how decay constraints may depend on the 
as-yet poorly-understood star-formation histories of the satellite galaxies. 

We use the prescriptions of Ref. [ |25| ] to estimate the minimum number of satellites with 
properties similar to the known population likely to be in the Milky Way halo based on the charac- 
teristics of the known satellites, taking into account the sky-coverage and completeness of existing 
searches for Milky Way satellites and the uncertainty in the Milky Way virial mass. We consider a 
decay model to be allowed if at least a few of our satellite-population realizations produce at least 
the minimum number of satellites we expect are in the Milky Way halo. 

Our constraints are summarized in Fig. ||], which shows the allowed and excluded regions 
of Vk — 1 parameter space assuming that the X particle behaves like cold dark matter. The pale 
red region to the right of the dashed curve and under the solid curve is excluded by the galaxy- 
cluster mass function and the mass-concentration relation. Local Group satellites provide weaker 
constraints in this region of parameter space. The pale red region to the left of the dashed line is that 
excluded by Milky Way satellites assuming that any subhalo of mass > 1O 7 M at its accretion time 
has not been merged with the main halo by dynamical friction; this is an extremely conservative 
constraint. The red region is the additional part of parameter space excluded by the Milky Way 
satellites if dynamical friction in the decaying-dark-matter cosmology is well-described by that in 
ACDM. The upper and lower solid lines correspond to the constraints if we assume that all subhalos 
are populated with stars, with the upper line corresponding to limits if dynamical friction operates 
as in ACDM, and the lower curve if dynamical friction is ineffective. 

The limits are not very sensitive to the star-formation prescription for satellite galaxies, but are 
sensitive to the evolution of subhalos and satellites after they are accreted onto larger halos. Better 
constraints on decay are only possible if this evolution is better understood. 



4. Discussion 



In this proceeding, we have described the constraints one may infer on decaying dark matter 
from different types of observations of dark-matter halos. The main result for decays is Fig. ||. 
The galaxy-cluster mass function and the mass-concentration relation are most constraining for 
Vk > 200 km s _1 , and the Milky Way satellites are most constraining for Vk < 200 km s , with 
X < 30 - 40 Gyr excluded for v k > 20 km s" 1 . 

There are several interesting ideas that have emerged as we have done this work. First, con- 
straints should improve greatly in the future with the advent of truly large and deep sky surveys 
[26]. Second, the cluster mass function and other probes of the growth of structure, traditionally 
used to constrain the normalization of the matter power spectrum and the nature of dark energy, 
can be used to constrain the nature of dark matter. However, the signature of decaying dark matter 
may be degenerate with that of modified gravity. This suggests that more work should be done to 
determine how well next-generation surveys can probe deviations of "A" and "CDM" simultane- 
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ously [17]. Third, the distribution of central densities of Local Group satellite galaxies should be a 



sensitive probe of the nature of dark matter-see Ref. [ |10| ] for further discussion. 
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